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Abstract 

We  propose  a  new  class  of  growth  and  interest  rate  processes  with  appealing  properties  for 
the  paper-and-pencil  theorist.  It  yields  simple  closed-form  solutions  for  stocks,  bonds,  and  perpe- 
tuities, and  can  accommodate  an  arbitrary  number  of  factors.  Expressions  are  linear  in  the  state 
variables,  such  as  the  interest  rate,  the  equity  premium  and  the  stock's  growth  rate.  Surprisingly, 
one  can  change  the  volatility  of  the  process,  or  force  the  interest  rate  to  be  positive,  without 
changing  bond  prices.  The  process  generaUzes  to  discrete-time  settings,  and  has  a  number  of 
economic  modeling  applications.  These  include  (i)  macroeconomic  situations  with  changing  trend 
growth  rates,  (ii)  asset  pricing  with  time-varying  risk  premia  or  time-varying  dividend  growth 
rates,  and  (iii)  yield  curve  analysis  that  allows  flexibihty  and  transparency.  (JEL:  G12,  G13) 

Keywords:  Long  term  risk.  Growth  rate  risk.  Perpetuity,  Modified  Gordon  growth  model.  Yield 
curve.  Interest  rate  processes,  Stochastic  Discount  Factor,  Quadratic  processes. 
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1     Introduction 

This  paper  proposes  a  new  growth  and  interest  rate  process  that  has  a  number  of  attractive  properties. 
Most  notably,  it  is  simple  and  tractable,  allowing  closed-form  solutions  for  the  prices  of  stocks,  bonds 
and  perpetuities.  It  is  also  flexible,  can  be  applied  to  an  arbitrary  number  of  factors  with  a  rich 
correlation  structure,  and  generalizes  well  to  discrete  time.  All  of  these  features  make  it  particularly 
accessible  to  the  paper-and-pencil  theorist. 

Quadratic  processes  have  the  drift  of  an  Ornstein-Ulenbenck  processes,  with  a  squared  term 
added.  This  squared  term  does  not  change  the  drift  much,  so  that  quadratic  processes  are  in  effect 
fairly  close  to  an  Ornstein-Uhlenbeck  process.  However,  adding  the  quadratic  terms  greatly  simplifies 
the  expressions  of  expected  values  of  stocks  and  bonds.  The  expressions  are  linear  in  the  factors. 
We  obtain  closed  forms  for  the  price  of  stocks  with  a  time-varying  growth  rate,  and  the  price  of 
perpetuities.  Such  closed  forms  for  perpetuities  is  not  not  available  with  the  available  processes,  such 
as  those  of  Ornstein-Uhlenbeck  /  Vasicek  (1977),  Cox,  IngersoU,  Ross  (1985),  the  Heath,  Jarrow, 
Morton  (1992),  or  models  in  the  affine  class  (Duihe  and  Kan  1996). 

Quadratic  processes  are  meant  to  be  a  practical  tool  for  several  topics  in  economics.  The  quadratic 
process  is  likely  to  be  useful  in:  (i)  macroeconomics,  for  models  in  which  GDP  growth  rate  is  time- 
varying;  (ii)  asset  pricing,  for  models  where  the  equity  premium  is  time-varying,  and  models  where 
the  trend  growth  rate  of  dividends  is  varying;  and  (iii)  fixed-income  analysis.  Also,  the  process  may 
be  useful  for  thinking  about  other  situations  where  the  uncertainty  in  the  discount  factor  matters, 
e.g.,  in  environmental  economics. 

Several  literatures  motivate  the  need  for  a  tool  such  as  the  quadratic  process.  Many  recent 
studies  investigates  the  importance  of  long-term  risk  for  asset  pricing  and  macroeconomics,  e.g., 
Bansal  and  Yaron  (2004),  Barro  (2006),  Croce,  Lettau  and  Ludivgson  (2006),  Gabaix  and  Laibson 
(2002),  Hansen  Heaton  and  Li  (2005),  Hansen  and  Scheinkman  (2005),  Julliard  and  Parker  (2004), 
Parker  (2001).  The  quadratic  process  offers  a  way  to  model  long-term  risk,  while  keeping  a  closed 
form  for  stock  prices. 

In  addition,  there  is  debate  about  the  existence  and  mechanism  of  the  time-varying  expected 
stock  market  returns,  e.g.,  Campbell  and  Shiller  (1988),  Cochrane  (2006),  and  many  others.  Because 
of  the  lack  of  closed  forms,  the  literature  relies  on  simulations  and  approximations.  The  quadratic 
process  offers  closed  forms  for  stocks  with  time-varying  equity  premium,  which  is  useful  for  thinking 
about  those  issues. 

Finally,  we  contribute  to  the  vast  literature  on  interest  rate  processes,  by  presenting  a  new, 
flexible  process.  It  is  as  flexible  as  the  afline  class  of  Duffle  and  Kan  (1996),  which  includes  the 
Vasicek  (1978)  and  the  Cox,  IngersoU,  Ross  (1985)  process  as  special  cases.  Section  4.3  develops  the 
link  between  the  quadratic  class  and  the  affine  class. 

This  paper  stipulates  a  process  for  finding  desirable  properties  for  the  pricing  kernel,  and  in  this 
follows  a  productive  literature  represented  by,  e.g.,  Campbell  and  Cochrane  (1999),  Cox,  IngersoU, 
Ross  (1985).,  and  particularly  Menzly,  Santos  and  Veronesi  (2004)  and  Pastor  and  Veronesi  (2005).^ 
The  last  two  papers  are  the  closest  to  this  paper,  as  they  contain  several  useful  closed  forms. 

The  core  insight,  that  the  quadratic  process  (72)  yields  particularly  tractable  formulas,  appears 
to  be  new  to  the  literature.'^ 


Those  papers  do  provide  an  economic  foundation  (via  external  habit  formation)  for  the  postulated  processes. 
The  Fisher- Wright  process  (e.g.,  Karlin  and  Taylor  1982)  does  contain  a  quadratic  term,  but  it  has  not  been  applied 


Section  2  presents  the  basic  one-factor  process.  Section  3  applies  the  one-factor  process  to  stocks. 
Section  4  presents  the  multifactor  version  of  the  quadratic  process,  which  is  useful  for  talking  about 
processes  with  fast  and  slowly  mean-reverting  components,  or  interest  rate  processes  with  short  and 
long  rates.  Section  5,  which  is  more  technical,  studies  the  range  of  admissible  initial  conditions. 
Section  6  presents  the  discrete-time  version  of  the  process.  Section  7  concludes. 

2     The  one-factor  quadratic  process 

2.1      Definition  and  basic  properties 

We  fix  a  probability  space  (fi,  T ^  P)  and  an  information  filtration  !Ft  satisfying  the  usual  technical  con- 
ditions (see,  for  example,  Karatzas  and  Shreve  (1991)).  We  start  from  a  simple  Ornstein-Uhlenbeck 
process,  but  add  a  quadratic  term: 

drt  =  -{P  -a)  {rt  -  a)dt+  {n  -  af  dt  +  dMt, 

with  a  </?,  for  rt  <  f  £  (a,/?)  where  rt  is  the  short-term  interest  rate. 

—  (/3  —  a)  {rt  —  a)  dt  is  the  standard  term  which  leads  to  mean  reversion  to  the  long  run  value 
a,  with  a  speed  (3  —  a.  The  second  term,  {rt  —  a)  ,  represents  a  small  departure  from  a  classical 
Ornstein-Uhlenbeck  process.  When  rt  is  not  too  far  from  a,  the  drift  is  little  changed  from  the 
standard  —  (/?  —  a)  (r  —  a)  dt^ .  However,  it  turns  out  that  the  correction  {rt  —  a)  dt  substantially 
improves  the  tractability  of  the  process. 

We  assume  that  the  martingale  {Mt,t  >  0}  has  properties  that  ensure  the  interest  rate  remains 
below  some  upper  bound  r  e  {a,P),  and  thus  does  not  explode.  One  example  is  dMt  =  o' {h)  dzt 
where  zt  is  a  Brownian  process  and  a  {r)  ^  k{r  —  r)'^,  for  r  in  a  left  neighborhood  of  f,  k  >  1/2 
and  A;  >  0.  Given  the  drift  is  negative  around  f,  that  will  ensure  that  f  is  a  natural  boundary,  and 
{Vf ,  rt  <  f}  almost  surely.''  We  also  have  the  option  to  set  a  lower  bound  r  to  the  interest  rate,  by 
adding  an  additional  restriction  on  dMt.  For  instance,  if  one  wishes  to  ensure  rt  >  r  for  some  r  <  a, 
we  could  have  dMt  —  o"  {rt,  t)  dzt,  with  at  ^  k'  {r  —  r)'^  ,  k'  >  1/2  for  r  in  a  right  neighborhood  of  0, 
and  k'  >  0. 

The  above  process  can  be  rewritten: 

drt  =  {rt  -  a)  {rt  -  P)  dt  +  dMt.  (1) 

Appendix  A  contains  a  formal  statement  and  justification  of  these  assumptions,  and  further  purely 
technical  assimiptions.  Given  these  assumptions,  and  an  initial  value  tq,  the  stochastic  differential 
equation  (1)  has  a  unique  strong  solution  for  rt,  defined  for  all  times  t  >  0. 

The  following  Theorem  gives  a  basic  result.  Its  proof  is  in  Appendix  B. 


to  the  pricing  bonds  or  stocks.  Also,  it  is  more  special  than  the  quadratic  class,  because  it  imposes  a  specific  functional 
form  on  the  variance.  Driessen,  Maenhout  and  Vilkov  (2005)  apply  the  Fisher-Wright  process  to  options. 

''in  many  applications,  it  is  plausible  that  \rt  —  q|  <  5%,  and  thus  the  quadratic  term  (rj  —  a)  is  less  than  0.25% 
per  year. 

■"See  Karlin  and  Taylor  (1981,  Lemma  15.6.3,  Table  15.6.2),  and  the  appendix  of  Ait-Sahalia  (1996)  for  a  concise 
treatment. 


Theorem  1    The  price  at  time  t  of  a  zero-coupon  bond  of  maturity  T,  Zt  (T)  =  Et   exp  (  —  J^       rsds 

is: 

Z.(T)  =  §^e-^+^e-'^^.  (2) 

p  —  a  p  —  a 


Expression  (2)  is  particularly  simple,  as  it  is  linear  in  the  current  rate  rt-  Remarkably,  the  bond 
price  does  not  depend  on  the  noise  dMt  and  thus  the  volatility  of  the  process.  ^  The  intuition  is  as 
follows.  Consider  the  case  where  the  process  is  deterministic,  i.e.,  dMt  =  0  for  all  t.  Theorem  1  still 
applies,  and  we  rewrite  Eq.  2  as:  Z  {T,r)  =  A  +  Br,  where  A  and  B  are  constants  that  depend  on  T 
but  not  on  r.  Suppose  that  at  time  t  the  interest  rate  is  rt,  and  apply  a  shock  to  the  interest  rate  at 
time  t  +  dt  so  that  rt+dt  =  ''"t+ut,  where  ut  is  a  mean-zero  random  variable.  After  t  +  dt,  the  process 
remains  deterministic,  and  evolves  according  to  (1)  with  dMt  =  0.  Given  an  initial  value  rt  +  ut  of 
the  interest  rate,  the  bond  price  at  t  +  dt  is  Z  (T,  rt  +  ut)  =  Z  {T,rt)  +  But,  because  Z  is  affine  in 
r.  Hence  the  value  at  i  is  £■  [Z  {T,  rt  +  Ut)]  =  Z  {T,  rt).  The  expected  value  of  the  impact  of  uj  is  0, 
because  the  bond  price  is  linear  in  uj.  A  zero  mean  addition  of  volatility  does  not  affect  the  bond 
price.  The  cnrx  is  that,  in  the  deterministic  case  {\ft,Mt  =  0),  the  bond  price  is  an  afhne  function 
of  the  interest  rate:    Z{T,rt)  =  A{T)  +  B  {T)rt.    By  contrast,  intuition  based  on  the  expression 


E, 


exp 


(-/; 


t+T 


„ds 


would  suggest  a  convex  function. 

'The  independence  of  bond  prices  from  volatility  greatly  simplifies  the  analysis.  In  particular, 
dMt  could  have  jumps,  which  model  a  decision  by  the  central  bank.  One  does  not  need  to  specify 
the  volatility  process  to  get  the  prices  of  bonds:  only  the  drift  part  is  necessary.  This  leaves  a  high 
margin  of  flexibility  to  cahbrate  volatility,  for  instance  on  interest  rate  derivatives,  a  topic  we  do  not 
pursue  here. 

Equally  surprisingly,  we  can  impose  a  lower  bound  without  affecting  bond  prices,  provided  the 
lower  bound  is  less  than  r<  a.  The  intuition  is  the  following.  We  can  make  r  a  lower  bound  for  the 
process  simply  by  setting  the  variance  to  be  0  below  r.  Since  variance  does  not  affect  bond  prices, 
the  new  lower  bound  does  not  affect  bond  prices. 

Existing  models,  such  as  Vasicek  and  Cox-Ingersoll-Ross,  are  able  to  generate  closed-form  so- 
lutions for  zero-coupon  bonds.  However,  the  quadratic  process  introduced  by  this  paper  can  also 
generate  closed-form  solutions  for  perpetuities. 


Proposition  1    The  price  of  a  perpetuity  is 


roo  .J, 


a  +  f3  -  rt  _  1       a  -  rt 
aP  a         a/3 


(3) 


Proof.  This  follows  directly  from  Eq.  2  and 


V 


roo 

/      Zt{T) 
Jo 


dT 


13  -rtl    ,   rt-  al       a  +  P  -  rt 


P  —  a  a       P  —  a  P 


ap 


We  now  analyze  the  dynamics  of  the  process  in  some  further  detail.  If  we  start  from  rt  =  a,  then 
Zt  (T)  =  e~°''^ ,  which  is  the  value  the  bond  would  have  if  rt  remained  pinned  at  a.   However,  rt  is 


'Hence  the  process  exhibits  "unspanned  volatility"  in  the  sense  of  ColUn-Dufresne  and  Goldstein  (2002). 


stochastic.  With  a  <  (3,  Eq.2  shows  that  Zt  (T)  ~  j^e~°''^  for  T  — >  oo.  Hence  a  is  the  long  run 
interest  rate,  a  =  —  limr-.cxj  Z'  (T)  /Z  (T).  Given  the  drift  process  (1),  E  [drt]  >  0  if  r  is  below  a, 
and  E  [drt]  <  0  if  r  is  above  a.  Hence,  r  is  mean-reverting  towards  a. 

P  has  two  roles.  First,  it  represents  an  upper  bound  to  the  interest  rate:  rt  <  f  <  p.  Second, 
P  —  a  is  the  speed  of  mean-reversion  of  the  process.  When  /3  is  high,  the  process  mean-reverts  faster. 

The  process  is  defined  only  for  rt  <r  <  j3.  If  r^  >  f,  the  process  may  explode  in  finite  time. 
Indeed,  when  the  process  is  deterministic  (M  =  0)  and  rj  >  f,  then  the  process  explodes  in  finite 
time,  i.e.,  there  is  ii  >  t  such  that  lims^fj  r^  =  oo.  The  assumption  r^  <  ^  is  needed  in  the  proof 
to  ensure  that  the  process  is  defined  at  times  in  \t,t  +  T].  This  is  guaranteed  if  rj  <  /3,  under  the 
conditions  of  Appendix  A. 

Finally,  we  require  the  leading  terms  to  be  r^  in  the  definition  (1)  of  the  quadratic  process,  not 
kr^  for  a  /c  7^  1. 

The  empirical  relevance  of  the  quadratic  term  is  unclear.  It  could  be,  for  instance,  that  the 
quadratic  term  exists  under  the  risk-neutral  probability,  though  not  in  the  physical  probability.  Also, 
of  course  all  processes  are  an  approximation  of  reahty,  and  in  many  cases,  for  instance  for  options, 
the  major  counterfactual  assumption  is  to  assume  Gaussian  innovations  rather  than  fat-tailed  ones 
(see  Gabaix  et  al.  2003,  2006).  Quadratic  processes  allow  for  jmnps  and  fat-tailed  innovations,  as 
the  shape  of  the  innovations  (the  dMt  term)  does  not  affect  the  prices. 

The  rest  of  this  section  considers  more  specialized  topics  in  the  one-factor  quadratic  process. 

2.2      Sum  of  two  quadratic  processes 

This  section  studies  the  sum  of  two  quadratic  processes.  This  apphcation  is  useful,  for  instance, 
when  both  the  equity  premium  and  the  dividend  growth  rate  are  time-varying. 

Proposition  2  (Sum  oj  two  independent  quadratic  processes)  Consider  two  one-factor  quadratic 
processes: 

drt     =     -cp{rt  -  rt)dt  + {rt -r^)'^  dt  +  dnt 
dRt     =     -^  {Rt  -  R»)  dt  +  {Rt  ~  R^f  dt  +  dNt 

where  n  and  N  are  martingales  with  uncorrelated  innovations,  d  (n,  N)^  =  0,  and  such  that  (n.  A'')  is 
a  martingale.   Then, 


Et 


exp 


t+T 


{rs  +  Rs)ds 


exp 


t+T 


r^ds 


Et 


exp 


t+T 


R.ds 


and  the  price  of  a  perpetuity  with  interest  rate  rt  +  Rt,  Pt  =  Et    /q°°  exp  (  —  ^^       (r^  +  Rs)ds\ 
is,  calling  S„  =  r^,  +  R^, 


Pt 


S. 


rt-r,       Rt-R.       {n  -  r,)  {Rt  -  R,)  {2S,  +  (f)  +  ^) 
S,  +  4>       5*  +  $         (5.  +  0)  (5.  +  *)  (5.  +  <^  +  $)   J 


(4) 
dT 

(5) 


Proof.  The  first  part  is  immediate  using  Ito's  lemma.  The  second  part  comes  from  the  calculation 
of: 


P,=  1^     e-(-+^-)^(l  +  (r,-r. 


-4>T  _  I 


l  +  {Rt-  R.) 


„-*T 


-1 


$ 


dT. 


3     Application  to  the  pricing  of  stocks  with  a  stochastic  growth  rate 


This  section  applies  the  quadratic  process  to  model  the  dividend  growth  rate  of  stocks.  We  start 
with  a  stochastic  trend  growth  rate  of  dividends,  continue  with  a  stochastic  equity  premium,  and 
finally  analyze  more  complex  examples. 

3.1      Time- varying  dividends 

The  stock  has  a  dividend  Dt,  which  grows  at  a  stochastic  rate  gt: 

^     =    gtdt  +  dNt  (6) 

dgt     =     ct>{g,-gt)dt-{gt-g,fdt  +  dMt  (7) 

(dMt,dNt)/dt    =    cgt  +  b  (8) 

9t     >    g>P  (9) 

where  P  <  a,  Nt  is  a.  martingale,  and  Mt  is  a  martingale  that  ensures  gt>g  for  some  g  >  p.  When 
6  =  c  =  0,  g*  is  the  long-run  growth  rate  of  the  stock. 

Innovations  to  dividends,  dNt,  can  be  correlated  with  innovations  to  the  growth  rate,  dMt-  The 
strength  of  that  correlation  can  vary  with  the  level  of  the  growth  rate  gt- 

The  discount  factor  is  constant  at  R.  The  stock  price  is: 


Vt  =  E 


e-^'^'-'^Dsds 
t 


We  start  with  the  case  b  =  c  —  0,  i.e.,  of  no  correlation  between  dividend  level  and  growth  rate. 

Proposition  3  (Modified  Gordon  growth  model,  with  time-varying  dividend  growth  rate)  When 
there  is  no  correlation  between  innovation  the  level  of  dividend  and  innovations  to  their  growth 
rate  ({dMt,dNt)  =  0),  the  price/dividend  ratio  of  the  stock  is: 

V^./A  =  ^  (l  +  ^^^^)  (10) 

R-  g*  \        R-  g^  +  4>J 

where  g^  is  the  long  term  growth  of  dividends,  and  gt  —  g*  is  the  deviation  of  the  current  growth  rate 
{gt)  from  the  trend  (g*). 

The  first  term,  R  —  g^,  is  the  P/D  ratio  from  the  Gordon  growth  model,  since  5,  is  the  long  term 
growth  rate  of  the  dividend.  The  second  term,  gt  —  g*,  is  the  deviation  of  the  current  growth  rate 
(gt)  from  the  trend  (5*).  It  increases  the  P/D  ratio,  by  an  duration  factor  1/  (/?  —  g,  +  (p)  that  is 
decreasing  in  the  discount  rate  R,  and  the  speed  of  mean  reversion  (p. 

We  now  consider  the  case  where  {dMt,  dNt)  /dt  =  cgt  +  b. 


Proposition  4   (Modified  Gordon  growth  model,  with  correlations  between  innovations  to  the  level 
and  growth  rate  of  dividends)  The  ■price/ dividend  ratio  of  the  stock  is: 


Vt/Dt 


R  -  g*  +  (p  +  c  +  {gt  -  g*) 
{R-g,){R-g,  +  (l>)-cR-b' 


(11) 


The  stock  price  is  increasing  in  the  correlation  between  the  level  of  dividend  and  its  trend  growth 


rate. 


3.2      Time- varying  equity  premium 

Consider  the  stochastic  discount  factor: 


Qt  =  exp  (  -ri  -   /    -^ds  -   /    X^dB, 


(12) 


The  price  at  time  t  of  any  asset  yielding  a  dividend  Dsds  at  time  s  is:   Vt  =  Et  [J^  QsDsds]  /Qt- 
We  seek  to  price  a  stock  with  dividend: 


Dt  =^exp(gt-   /    ^ds  +   /   a^dBs )  Dq. 


(13) 


Assume  that  the  risk  premium,  nt  =  XtO't,  follows  a  one-factor  quadratic  process: 

dnt  =  -(p  {-Kt  -  TT,)  +  (ttj  -  TT,)^  dt  +  dNt 

where  Nt  is  a  martingale  that  ensures  that  the  process  is  well-defined  (tt^  <  tF  <  tt*  +  0)-  Assume 
that  processes  TV^  and  Bt  are  independent.  Then,  the  price  of  a  stock  is,  calling  J^/^  the  filtration 
associated  with  Nt, 


QtDtdt 


Uo 


/Qo 


/      exp  Hg  -  r)  i  -    / 


Vo  =  Eo 
=  Eo 
=     Eo 

=    ^■o 
so  that,  applying  Proposition  1  to  rt  =  —g  +  r  +  nt,  we  get: 

1  /.  TTt  -  TT* 


oo  /  rt    \2 

expUg-r)t-  ^ 


*  \2    ,      2  ft 

'^    'ds+        ias-Xs)dBs]dt 

rt  /_  ^    ^2 

ds  + 

---      Jo 


Xi  +  a',,     ,     r{as-Xs)   ^^U,|^N 


/      expl{g  —  r)t—   /    Xsasdsjdt    =  Eq     /      exp{{g  —  r  —  Trt)t)dt 


Vt/Dt 


r  +  7T»  -  g 


1 


r  +  7rt  -  g  +  ' 


(14) 


This  is  a  Gordon  growth  formula,  with  time- varying  risk-premium  ttj.   When  the  risk  premium 
TTf  is  equal  to  its  central  value,  tt*,  we  get  the  traditional  formula,  Vt/Dt  —  1/  (r  -|-  tt,  —  g).  When  ttj 


is  different  from  tt,,  we  get  an  adjustment,  equal  to  ttj  —  tt,  times  the  duration  of  that  term,  which  is 
the  effective  rate  of  discounting,  r  +  tt*  —  5,  plus  the  speed  0  of  mean-reversion  of  the  risk-premium. 

3.3      Time- varying  equity  premium  and  dividend  growth  rate 

We  now  merge  the  two  previous  examples,  to  incorporate  both  a  time- varying  equity  premium  and  a 
time-varying  dividend  growth  rate.  The  stochastic  discount  factor  and  dividend  are  given  as  follows: 


ft  )2  rt 

Qt  =  exp  {-rt-   /    -^ds  -   /    XsdBs 


(15) 


(/' 


gsds  -  -^ds  +  asdBs  )  Dq, 


Dt  =  exp 
where  gt  follows  the  quadratic  process 

dgt  =  -$  (gt  -  g*)  dt  -  {gt  -  g*f  dt  +  dNt. 
The  risk  premium,  ttj  =  XtCti  follows  a  one-factor  quadratic  process: 


d-Kf 


-(j)  {iXt  -  TT*)  dt  +  (TTt  -  TT,)     dt  +  dnt 


where  nt,  Nt  are  martingales  satisfying  the  familiar  quadratic  conditions.  Assume  that  the  processes 
nt,Nt  and  Bt  are  independent.  Then,  the  price  of  a  stock,  Vt  =  Eq  [J^  QtDtdt]  /Qq,  is,  by  the 
reasoning  of  the  previous  section. 


Vt  =  Et 


exp 


L"'5=f 


{r  +  Tru-  9u)  du  1  ds 


u=t 


Define  5,  =  r  -|-  tt,  —  5*.  Formula  (5)  gives: 


V.ID.  =  ^ 


1  _  ^i  ""'^*     ,    3t 


*  _  (TTt  -  ttQ  {gt  -  g,)  (25.  +  cj)  +  ^) 

5*  +  $        {S.  +  (j>)  (5*  +  ^){S,  +  4>  +  $) 


(16) 


The  central  value  is  again  the  Gordon  formula,  Vt/Dt  =  1/  (1  +  tt,  —  g^).  It  is  modified  by  the 
current  level  of  the  equity  premium,  and  the  growth  rate  of  the  stock.  A  stock  with  a  currently  high 
growth  rate  gt  exhibits  a  higher  price-dividend  ratio,  and  this  is  amplified  when  the  equity  premium 
is  low,  as  shown  by  the  term  (ttj  —  T^*){gt  —  9*)- 

In^the  above  process,  ^("ffleaTFreverts,  so  the^rice/dividend  ratio  mean-reverts  in  (11).  There 
is  much  recent  literature  on  the  predictabiUty  of  returns,  in  particular  the  predictabihty  of  the  P/D 
ratio  (Campbell  and  Shiller  1988;  Cochrane  2006).  The  above  setup  may  serve  as  a  useful  benchmark 
to  evaluate  the  debate. 


^Menzly,  Santos  and  Veronesi  (2004,  Eq.  20)  have  a  similar  expression,  even  though  the  processes  are  quite  different 
in  their  details. 


3.4      Correlation  between  the  dividend  and  the  discount  factor 


We  next  study  an  example  where  the  innovations  to  the  dividend  are  correlated  with  the  discount 
factor.  We  suppose: 

-—^  =  gdt-V  dNt 

with  7  is  a  fixed  number,  and  Nt  a  martingale.  We  assume  a  quadratic  process:  drt  =  [rt  —  A)  (r^  —  /x)  dt+ 
dMt,  and 

'  dDt    ,    \  ,  , 
-—-,drt)  /dt  =  an  +  b 

for  two  constants  a  and  b.  Dividend  growth  may  be  correlated  with  the  interest  rate. 

Proposition  5    The  price  at  time  t  of  a  claim  at  date  s  >t  is: 

E  \e-  f'  -'^"D  J  =  f  e-^'(-')  4^  +  e-^'(^-')  "-^)  D,  (17) 

L  J\  ^  —  X  ji   -  \  J 


where  A'  and  ^'  satisfy 


X'  +  ^'     =     X  +  /J.  —  g  —  a 
AV     =     {X-g)iix-g)  +  b 


(18) 
(19) 


i.e.,  are  the  solution  of 


e-lX  +  ii-g-a](  +  iX-g){i.L-g)  +  b  =  0. 
Proposition  6    The  price  at  time  t  of  a  dividend  claim  yield  Ds  at  every  s  >  t  is: 


E 


,-i:ruduj^    , 


X  +  ii-g-a-rt 
{X-g)ifi-g)  +  b 


Dt. 


LJs=t 
Proof.  By  direct  integration  of  (17).    ■ 

3.5      A  quantitative  illustration 

When  gt  is  at  its  long  run  value,  g^,  Eq.  10  gives  the  simple  price  dividend  ratio  1/  {R  —  gt).  However 
the  P/D  ratio  also  varies  with  the  growth  rate.  To  assess  the  quantitative  importance  of  this  effect, 
consider  the  following  calibration  using  aimualized  values:  R  =  8%,  which  reflects  a  risk-free  rate 
of  2%  and  an  equity  premium  of  6%;  g^  =  3%,  the  long-run  growth  rate,  which  is  proportional  to 
typical  values  of  long-run  GDP  growth.  This  yields  a  central  P/D  ratio  of:  1/  {R  —  g,)  =  20.  We 
also  take  </>  =  5%,  which  gives  a  half-life  of  mean  reversion  of  \ii2/4>  =  14  years.  This  gives 


Vt/Dt  = 


1 


+ 


9t 


R-g.       {R  -  g.)  {R  -  g,  +  4>) 


=  20-1-  200  {gt  -  a) . 


With  a  volatility  ag^  —  2%,  we  get  a{V/D)  —  200(7g  =  4.   The  P/D  ratio  has  an  annual  standard 
deviation  of  4,  and  the  volatility  of  returns  is,  at  gt  =  g* 

-.„(v./oo  =  ^3^  =  10-.  =  20%. 

3.6      Application:  Long  term  gro^vth  risk 

Consider  dCt/Ct  =  gt,  with  gt  stochastic.  Thus  Ct+r/Ct  =  exp  (J^       gsds).  The  price  of  a  stock  in 
such  a  model  is: 

Vt/Ct  =  J      ^^'^i-^]        dT=j      exp(-p5-(7-l)y        gsdsUT. 

If  we  define: 

n  =  P+(7-  l)5i, 

we  obtain  Vt/Ct  =  /q°°  e^Jt      '^"'^^ds,  the  perpetuity  price. 

We  study  the  process  in  two  representative  cases.  For  a  quadratic  process,  drt  =  [rt  —  A)  (rj  —  /z)  (it+ 
(Ttdzt,  with  A  =  /O  +  75*  <  /i  we  have: 


A  \  /U 

Therefore,  the  volatility  of  the  P/E  ratio  is,  at  the  central  value  rt  =  A: 

0"r|r=A 


'^dlVi/Wt 


/^ 


We  calibrate  using  7  =  3,  5  =  2%,/)  =  1%,  which  yields  A  =  7%.  We  choose  a  central  P/E 
ratio  of  1/A,  which  is  14.3.  Around  A,  the  process  behaves  like  a  Ornstein-Uhlenbeck  process, 
dr  =  —  (/J  —  A)  (r  —  A)  dt  +  cr\dzt.  We  select  jj.  =  15%  and  a  speed  of  mean-reversion  jj,  —  X  =  8%, 
i.e.,  a  time  scale  of  1/  {jj,  —  X)  =  12.5  years. 

We  pick  (Jr  =  0.5%.  This  yields  a  standard  deviation  of  growth  rate  of  approximately:  ar  {2{fj,  —  X))     ' 
ar  (2  ■  8%)~  '     =  (Tj./0.4.     Taking  1.5%  as  the  standard  deviation  of  the  growth  rate,  we  have: 
ar  =  0.4  ■  1.5%  =  0.6%.  With  a  growth  rate  volatility  of  ag  =  0.5%,  we  obtain  interest  rate  volatihty 
of  ar  =  "fCTg  =  1.5%.  Choosing  fj.  =  15%,  this  yields  a  stock  price  volatility  aw  =  y^  —  10%. 


4     The  multifactor  quadratic  process 

Several  factors  are  needed  to  capture  the  dynamics  of  bonds  (Litterman  and  Scheinkman  1991)  and 
stocks  (Fama  and  French  1996).  Accordingly,  we  study  the  multifactor  version  of  the  quadratic 
process. 
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4.1  Definition  and  basic  properties 

A  quadratic  process  X^  S  M"  has  the  form 

dXt    =    {b  -  AXt)  dt  +  {n  -  r^)  (Xt  -  X,)  dt  +  dMt  (20) 

n     =    /3'(Xt-X.)  +  r,  (21) 

with:  yl  is  a  n  X  n  matrix  with  positive  eigenvalues  {5i)-_-^  „.  b,P,  G  ]R",X,  =  A~^b  e  lR",r*  £  R. 
Mt  6  M"  is  a  martingale,  so  E[dMt]  =  0,  but  its  component  dMu,  dMjt  can  be  correlated.  As  in 
the  one-factor  process,  the  volatility  of  dMt  must  go  to  zero  in  some  limit  regions  for  the  process  to 
be  well-defined.  We  defer  this  more  technical  issue  until  later. 

We  calculate  the  price  of  a  zero  coupon,  Zt  (T)  =  Et   exp  ( —  J^       rgds  j   . 

Theorem  2  (Price  of  a  bond  for  a  n— factor  process)  The  price  at  time  t  of  a  zero  coupon  bond  of 
maturity  T,  given  the  state  Xt,  is: 

p-AT  _  1 
Zt  (T)  =  e-^*^  +  e-'-' V J {Xt  -  X,) .  (22) 

We  use  a  function  /  (A)  for  a  matrix.  The  notation  is  standard.  If  /  is  analytic  and  admits  the  ex- 
pansion f  (x)  =  Yl  fn^^,  then  /  (A)  =  J2  InA'^.  If  A  is  in  diagonal  form,  A  =  QDiag  {5i, ...,  5n)  ^~^, 
then  /  (A)  =  QDiag  (/  (Ji) , ...,  /  (5„))  fi-^. 

Proposition  7  (Price  of  a  perpetuity  for  a  n— factor  process)  The  price  of  a  perpetuity,  Pt  = 
L    Zt  (T)  dT,  is,  with  In,  the  identity  matrix  of  dimension  n: 

Pt  =  --  -/3'  {A  +  rJn)-'  {Xt  -  X,) .  (23) 

4.2  Some  exEimples 

Example  1.     For  n  =  1,  we  obtain  the  expressions  of  the  one- factor  quadratic  process. 

Example  2.  Dividend  growth  rate  as  a  sum  of  mean-reverting  processes  (e.g.,  a  slow 
and  a  fast  process) 


gt    =   ff*  +  X]  -^^^ 


dXit    =    -(piXudt  -  {gt  -  g*)Xitdt  +  dMu 

The  growth  rate  gt  is  a  steady  state  value  g^,  plus  the  sum  of  mean-reverting  processes  Xu-   Each 
Xit  mean- reverts  with  speed  Ai,  and  also  has  the  quadratic  perturbation  {gt  —  g^)  X^dt. 

We  apply  the  Theorem  2,  with  P'  =  (1,...,1),  A  =  Dia^  ((;/>i, ...,  0„).    The  expected  value  of  a 
dividend  T  periods  ahead  is: 

"    1-e-*--^ 


Et  [Dt^T]  /Dt  =  e^-^  +  e^'^  ^        ^         X 


1=1 
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and,  for  the  price-dividend  ratio  of  a  stock: 

^'/^^  =  ^  f  1  +  E  .     f\.]  ■  (24) 

R-9*  \        ~{R-  9*+9i) 

Each  component  Xit  perturbs  the  baseline  Gordon  expression  1/  {R  —  g*).  The  perturbation  is  Xit, 
times  the  duration  of  Xi,  discounted  at  rate  R  —  g^,  which  is  the  term  1/  (i?  —  g*  +  (^J. 


Example  3.   A  short  rate  and  a  long  rate 


dr, 


cf>{Lt-rt)  +  {rt-r,f 


dt  +  dM} 


2 


dLt     =     [c^{r^-Lt)  +  {rt-r^){Lt-r^)]dt  +  dM\ 

Tt  is  the  short  term  rate.  Lt  is  the  long  term  rate.  Because  of  the  first  equation,  rj  is  drawn  to  Lt- 
The  second  equation  expresses  that  Lt  goes  to  a  steady-state  value  r,.  The  process  is  enriched  by 
the  terms  {rt  —  r,)  and  (rj  —  r*)  {Lt  —  r»),  which  make  the  process  quadratic,  hence  making  bond 
prices  tractable.  In  practice,  it  is  reasonable  to  expect  \rt  —  r^\  and  \Lt  —  r,|  both  less  than  5%,  so 
that  the  magnitude  of  the  quadratic  terms  is  less  than  5%^  =  0.25%  per  year. 

We  apply  the  Theorem  2,  with  Xi'=(^'j,/3=(M,A=(^  "J'  r^  ^  (  L  I'  '^'^ 
get,  with  r't  =  rt-  r*,  L't  =  Lt  -  r,, 


(j)        ip  J  ip  (j)  —  ip  \(j)       (p  —  ip 


Pt     ^     ~  '''  "^^'^ 


r,       r.  (r,  +  (j>)       r,  (r,  +  4>)  (r,  -|-  ip) ' 

Example  4.  rt  having  a  time-varying  trend  In  the  post-Volcker  era,  interest  rates  tended 
to  have  predictable  trends  of  increase  or  increase,  which  may  be  captvued  by  the  following  trend 
growth  rate  Sj  of  the  interest  rate: 

drt     =    stdt  +  in  -  r^f  dt  +  dMu 

dst    =     [-XfJ.{rt- r^,)  -  {X  + iJ.)st]dt+ {rt-r^)stdt  +  dM2t 


with  A,  /U  >  0,  and  A  -|-  ^  >  0.  Economically,  Sj  is  the  predicted  trend  in  interest  rates,  as  per  the 
first  expression.  St  mean-reverts  for  two  reasons:  first,  because  of  the  —X/j,  {rt  —  r,)  term  {st  becomes 
negative  if  interest  rates  are  too  high);  second,  because  of  the  —  {X  +  ^)  st  term. 
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We  apply  Theorem  2,  with  Xi  =  {n,  st),  /?'  =  (1,  0),A={^      ^  ^      ) .  We  get: 


Z{T)     =     e-^'^ 
Pt    = 


V  /        A  (/i  —  A)  '       fj.[fj.  —  \) 

1        St  +  (r,  +  A  +  ^)  (rt  -  r») 


r,  (r*  +  A)  (r.  +  ;u) 

Those  examples  show  it  is  quite  easy  to  get  closed  forms  with  easy  to  interpret  processes. 

4.3      Comparison  \vith  existing  processes 

Comparison  with  the  affine  class     The  afhne  class  (Duffie  and  Kan  1996;  Duffie,  Pan  and 
Singleton  2000;  Dufhe  2002)  comprises  processes  of  the  type: 

dXt     =     {b-  AXt)  dt  +  wtdzt 
wtw[     =    a^{H[Xt  +  Ho) 

with  Tt  =  u+  0'  {Xt  -  X,).  b,Xt  e  M",  A  e  R"^",  {Ho,Hi)  e  M"  x  ]R">^",cr  e  E.  We  note,  as 
before,  X,  =  A~^b,  assumed  to  exist. 

Under  mild  technical  conditions,  bond  prices  have  the  expression: 

Z^^  (T)  =  exp  {-r.T  +  T  {T)'  {Xt  -  X,)  +  a^a  (T)) 

with  r, where  a  (T)  and  T  (T)  satisfy  a  typically  compUcated  ordinary  differential  equation.  The 
situation  is  simpler  if  Hi  =  0.  In  that  case,  T  (T)  =  7  (T),  with  7  (T)'  =  P'  [e'^'^  -  l)  M-  Then: 

Z^^  {T)  =  exp  {-r.T  +  7  (T)'  {Xt  -  X,)  +  a^a  {T))  .  (25) 

The  express  can  be  contrasted  with  the  expression  for  the  quadratic  process  (20): 

Z^^{T)  =  e-^'^{l+'r{TyiXt-X.)).  (26) 

If  7  (T)  Xt  is  small,  the  two  expressions  are  the  same,  up  to  terms  of  second  order  in  7  (T)  Xt,  and 
second  order  in  a.  Hence,  a  quadratic  process  is  a  good  approximation  if  the  underlying  process  is 
in  fact  affine,  and  vice-versa.  In  most  cases,  the  two  values  are  likely  to  be  close,  so  that  existing 
estimates  of  parameters  in  the  affine  class  can  be  used  to  calibrate  quadratic  processes.^ 

Comparison  with  Heath,  Jarrow  and  Morton     In  the  Heath,  Jarrow  and  Morton  model, 
forward  rates  are  factors.    In  the  quadratic  process,  bond  prices  are  factors.    The  HJM  does  not 


'That  equivalence  gives  a  useful  way  to  calculate  easily  functionals  of  quadratic  processes,  that  can  be  expressed 
as  a  linear  combination  of  bonds.  One  first  works  with  the  affine  process,  setting  volatility  to  0,  doing  a  first  order 
Taylor  expansion  of  terms  in  {Xt  —  X.).  One  gets  an  expression;  P  =  a  +  b{Xt  —  X,)  +  o{Xt  —  X,)  -'r  o\u  ), 
for  some  constant  a,  6.  Then,  one  knows  that  for  the  corresponding  quadratic  process,  the  value  of  the  asset  is; 
P^^  =a  +  h{Xt-X.),  exactly. 
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yield  essentially  any  closed  forms,  but  has  proven  very  useful  to  price  options.   Applications  of  the 
quadratic  process  to  options  are  not  investigated  in  this  paper. 

4.4      Additional  remarks 

The  following  relation  highlights  why  the  quadratic  process  yields  linear  expressions. 

Lemma  1  Suppose  that  process  yt  G  M.'' satisfies  Et  [dyt]  /dt  =  Cytdt  +  rtytdt,  for  some  q  ^  q  matrix 
C.   Then,  VT  >  0, 

'        -    ^'■'^'^  T    ds 


Et 


yt+T 


e'^^yt. 


(27) 


Proof.  Call  Ls  =  Et 


-ffr^ds 


ys 


dLs  =  Et 


^-ffr,ds  ^_rsysdS  +  CysdS  +  rsysdS) 


CLs 


As  Lt  =  yt,  we  get  Ls  ^  e^^  ^^'^yt.    m 

As  a  corollary,  consider  which  would  have  a  constant  price.  We  can  easily  see,  using  the  above 
Lemma,  that  an  asset  yielding  a  dividend  Dg  =  r*  (l  —  P'A~^Xs)  at  each  date  s,  has  a  constant 
price  of  1. 

We  close  with  the  following  Conjecture.  If  true,  the  quadratic  term  is  the  only  way  to  obtain 
linear  bond  prices. 

Conjecture  1  Suppose  that  Xt  is  a  Markov  process  with  values  in  E"  such  that  the  interest  rate  is 
a  function  r  —  p  {Xt),  and  bond  prices  are  affine  in  Xt,  i.e.,  there  are  deterministic  functions  a  (T) 
and  r  (T)  with  values  in  M  and  R"  respectively,  such  that,  for  all  X  in  the  domain  of  definition, 

\/T  >  0,E 
process. 


exp  (  —  Jj       p{Xs)ds]  \  Xt  —  X     =  a(T)  +r{T)'X.     Then,  Xt  follows  a  quadratic 


5     Conditions  to  keep  the  process  well-defined 

5.1      The  practical  bottom  line 

With  one  factor,  the  process  is  well-defined  if  it  stays  within  r  <  f,  with  f  <  fx.  For  an  n— factor 
process,  the  following  algorithm  gives  the  conditions  for  the  initial  state  vector  to  yield  a  well-defined 
process. 

First,  one  diagonalizes  the  matrix  A,  i.e.,  sets  A  =  QAQt^,  with  A  =  Diag{5\,  ...,5„),  with  the 
root  not  necessarily  ordered.  The  eigenvector  corresponding  to  eigenvalue  5j  is  (Oij)^_j    ^. 

Then,  one  sets:^  Q  =  Diag  ({n'l3)j  /6j]  ■  fi-^ 

The  process  is  well-defined,  if  the  initial  value  of  Xt  satisfies  the  n  inequalities: 


Vi=  L..n,^ 


mm{5i,5j] 


Qjk  {Xk  -  X^k)  <  1- 


(28) 


This  rescales  the  eigenvectors,  so  that  in  the  diagonal  representation,  r  =  5' X . 
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which  we  call  define  being  in  the  "Root"  region. 

In  matrix  form,  one  defines  Tij  =  ( '"'"^^''  ''  j ,  and  the  criterion  is:  FQ  {Xt  —  X«)  <  (1, ...,  1),  in 
the  sense  that  the  inequality  holds  for  each  of  the  n  coordinates. 

If  the  initial  value  of  Xt  satisfies  (28),  and  increments  are  continuous,  then  all  future  Xgyt  also 
satisfy  (28).^ 

The  above  conditions  are  simply  sufficient,  but  one  conjectures  that  they  are  actually  optimal, 
in  the  domain  of  conditions  requiring  only  n  linear  inequalities. 

Some  examples  For  n  =  1,  the  condition  is  just  rt  <  fj,,  i.e.  rt  —  r,  <  (f>.  The  equivalent 
condition  for  the  stock  with  time- varying  growth  rate  is:  gt  —  g*  >  —(p- 

With  n  factors,  if  the  process  is  diagonal,  i.e.,  A  =  Diag  {6i,  ...,5n),  with  r  —  Y2Pi^i^  ''hen  ^ 

is  the  identity  matrix,  Q  =  Diag  (^Pj/5j),  and  the  condition  is  :  Vi,  ^"_j  """^"  ^  -^Xj  <  1,  i.e., 

w '   \~^^  ^jj        ^  1 

V^2^j=l  max((5.,(5j)   '^  ^■ 

When  n  —  1,  with  A  =  Diag  (61,62)  and  r  =  J2i=i  ^i^ii  the  root  region  is: 

X  I  Zi  +  X2  <  1  and  -iXi  +  X2  <  1 
62 

In  Example  2  of  section  4.2  with  a  dividend  growth  rate  gt  =  g*  +  ^1^=1  -^it,  ,  one  applies  the 
criterion  to  rt  =  i?  —  gt,  with  Pi  —  —1.  Therefore,  the  condition  is: 

Vz  =  l...n,V ^^>-l. 

^-r'  max(0i,0j) 

In  Example  3  of  section  4.2,  with  the  short  rate  and  the  long  rate,  we  get  5'  =  {(f>,ifj),  fl  = 
J,  )i  Q  =  Diag  {\/(j),(f>/'ip)Q,~^ .  When  cp  >  ip  (the  short  rate  moves  faster  than  the  long 

rate,  so  to  speak),  T  =  (    -,        (       I ,  and  the  conditions  become: 

— ■ —  <  1  and  — ; \ ; —  <  1. 

0  (p  ip 

Intuitively,  r  and  L  cannot  deviate  from  their  central  value,  r, ,  by  a  value  order  of  magnitude  the 
speed  of  their  mean-reversion. 
__  .    The  rest  of  this  section  justifies  these  claims.  It  may  be  skipped  in  a  first  reading.        _     _        __ 

5.2      Notation  and  motivations 

We  consider  a  diagonal  process,  i.e.,  A  =  Diag{5i,  ...,(5„),  with  6\  <  ...  <  5n,  and  n  >  1.  In  a  dense 
open  subset  of  cases,  we  can  rely  on  this  case  after  diagonalizing  A.  We  define  (5o  =  0  and  5n+i  =  6n- 
We  rescale  X  if  needed,  so  that  P  =  6,  and  so  that  r  =  ^  6iXi.  We  reason  first  in  the  deterministic 
version  of  the  process. 


One  can  show  that  another,  strictly  more  stringent,  sufficient  condition,  is  Vi  =  l...n,  X^"=i  ^Si<Sj  Qjk  {^k  —  X.k)  < 
1. 
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We  use  n  tests  vectors,  F^^^  ...,  r("\  with  Vi  =  1...,  n,  F'*^  e  M".  We  study  regions  of  the  type: 

wr  =  {x  G  M"  I  Vz  =  l...n,  F^')  ■  X  <  l| .  (29) 

We  want  to  make  sure  that,  if  the  process  hits  the  boundary  of  to,  noted  duj,  then  it  goes  back  inside 
w.  To  express  this,  call  m{X)  =  E  [dXt/dt  \  Xt  =  X]  =  {r  -  A)  X  the  drift  of  Xt  at  X.  If  X  £  du, 
then  call  j  the  coordinate  such  that  F'-^^  ■  X  =  1;  then,  we  want  to  ensure  F'^^  ■  m[X)  <  0.  We 
calculate: 

f(^)  •  to  (X)     =    f(^)'  {r-A)X  =  rV^'^'X  -  T^'^'AX 

=    r-1-  (aT^'A' X=  (S-AT^'A-X. 


So,  we  require,  for  all  i  >  1:  {Vj  ^  i,  F^^')  •  X  <  1  and  F^^)  ■  X  =  l}  =^  {S  -  AF^'))  ■  X  <  0. 

We  impose  this  condition  only  for  i  >  1,  because  in  practice,  a  boundary  of  choice  is  F'^'  = 
(1, ...,  1)',  which  requires  a  special  treatment.  Call  St  =  F'^^  ■  Xt  =  J2j  -^jt^  then  dSt/dt  =  rj  {St  —  1). 
So,  St  cannot  cross  1,  and  if  St^^  <  1,  at  an  initial  date  St^  <  1,  then  for  all  i  >  to,  S'f  <  1.  So,  the 
region  {X  e  R"  |  F'^^  ■  X  <  l}  is  automatically  non-repelling. 

The  following  Lemma  completes  the  reasoning,  by  providing  a  sufficient  condition  for  the  region 
iOr  to  be  non-repelling. 

Lemma  2  (Sufficient  condition  to  get  a  non-repelling  region)  Suppose  T'^^'  =  (1, ...,  1)',  and  that  for 
each  i  =  2...n,  there  is  a  ai  >  0  and  a  non-negative  vector  ip^^'  G  M",  with  ^  ■  ipj  =  1;  Vi  <  1; 
such  that: 

ai5+{l-  aiA)  F^')  =  ^  i>fr^^\  (30) 

j 
Then,  the  process  does  not  escape  from  region  w. 

Proof.  Indeed,  then 


a,  [5  -  AF«)    X    =     [  -F(^)  +  J2  ^f  r(^)      •  X 


5.3     The  Root  region 


-.(i) 


We  consider  the  region  wr  defined  by  F^.     =  5i/s,j/5i,  where  i  Aj  =  min(i,j).    Because  it  involves 


J 


the  eigenvalues  of  A,  we  call  it  the  "Root  region".  Heuristic  arguments  lead  to  the  conjecture  that 
the  Root  region  is  optimal,  in  the  sense  that  it  may  be  the  largest  non-repelling  region  that  can  be 
defined  by  the  intersection  of  no  more  than  n  inequalities  expressed  as  linear  combinations  of  the 
factors. 
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Theorem  3    The  Root  region  wp,  defined  by  the  boundary  vectors  Tj    =  5i/^j/5i  and 

0Jr=  {x  eW\yi=  l...n,r(^'-X  <  l| , 
is  non-repelling. 

5.4      Conditions  on  volatility 

In  dimension  1,  we  require  that  var  [dMt)  =  0  near  the  boundary.  It  can  be  0  in  a  region  [//  —  e,  /x], 
for  some  £:  >  0,  or  o"  (r)  ~  A;  (f  —  r)°  with  a  >  1/2  or  a  =  1/2  and  k  large  enough. 

In  dimension  n,  the  analog  is  that  when  r('''X(  =  1,  (while  X  €  aJr),  then  var  (r''^'  •  dMt) 
should  be  0  near  the  boundary.  The  simplest  condition  is  that,  for  each  i  =  l...n,  there  is  an  open 
neighborhood  i^(')of  wfH  [X  \  Mi,T^^^  ■  X  =  l],  such  that 

Xt  G  v^^'  =^  var  (r'^''  ■  dMt)  =  0.  In  other  terms,  the  noise  along  the  projection  on  T^^'  becomes 
0. 

A  a  simple  way  to  ensure  that  the  process  is  well-defined  is  the  following.  For  an  e  >  0,  in 
practice  very  small,  we  define  the  killing  function 


Ke  (X)  =  1  I  max  T^^  •  X  <  1  -  ej 


(31) 


where,  as  always,  1  {A)  is  equal  to  1  if  A  is  true,  or  0  otherwise. 

We  start  from  an  arbitrary  continuous  martingale  process,  dMt,  and  then  define  a  new  continuous 
martingale: 

dMt  =  K,{Xt)dMt.  (32) 

In  other  terms,  aU  variances  become  0  if  X  is  within  a  sraaU  neighborhood  of  the  limit  of  the  root 
region  wp-  Otherwise,  the  process  is  not  modified.  This  "killing"  procedure  might  be  the  most 
convenient  to  use  in  practice. 

5.5      Condition  with  diagonalizable  matrices 

In  the  general  case,  given  a  diagonalizable  matrix  A,  we  diagonalize  it,  A  =  nAn~\  with  A  = 
Diag  {5i,  ...,(5n)-  We  set  Q  =  Diag  ((J7'/3)j  /5i)-Q.''^,  and  define  the  new  state  vector  to  be:  Yt  =  QXt- 
Then: 

5'Yt    =    5' ■Diag{{n'(3)J5i)n-^Xt  =  Diag{{a'f3)^)n-'^Xt 

=    l3'n-n-^Xt  =  l3'Xt  =  rt. 

As  rt  =  5'Yt,  Yt  satisfied  the  conditions  of  this  section  for  the  tests,  and  the  Root  region  is:Vi  =  l...,n, 
r(')  •  y  <  1.   As  r(')  •  Yt  =  r'-'^'QXt  =  Ylj  ^QjkXk,  the  condition  can  be  re-expressed  in  a  way 

that  does  not  depend  on  the  prior  ordering  of  the  roots:  Vi,  ^  '^  ^  QjkXk  <  1,  which  gives  the 
condition  announced  in  (28). 
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6     The  process  in  discrete  time 

This  section  studies  the  discrete-time  version  of  the  quadratic  process,  which  is  useful  for  several 
reasons.  In  many  applications,  discrete-time  is  simpler  than  continuous  time.  The  technical  issues 
are  clarified  by  discrete-time  examples. 

The  key  object  is  the  stochastic  discoimt  factor  mt  6  K!^.  The  price  at  time  t  of  a  claim  yielding 
a  stochastic  dividend  Ds  at  date  S  >t  is:  Pt  =  E  [mt+i...msDs]-  The  price  of  a  zero  coupon  bond 
of  maturity  T  is: 

Zt  (T)  =  Et  [mt+i...mt+T]  ■  (33) 

For  instance,  the  gross  short  term  interest  from  t  to  i  -|- 1  is  l/Et  [mt+i\-  In  terms  of  the  continuous 
time  process,  mt  =  exp  f  —  j^_.^  Vudu  1 . 

6.1      Discrete  time,  one  factor 

The  most  obvious  generalization  of  the  process,  Emt+\  quadratic  in  tti^,  does  not  work.  Instead,  the 
generalization  we  propose  is: 

Etmt+i  =  -  +  \-  -^  (34) 

a       0       abrrit 

with  0  <  a  <  6.  We  assume  that  the  process  is  such  that  for  all  times,  mt  >  1/b}^ 

The  expected  value  of  m,t+i  is  increasing  in  mt,  so  the  process  has  a  positively  autocorrelation.  On 

the  other  hand,  if  there  is  no  noise,  then  m,t  goes  to  the  attractive  root,  1/a.  Indeed,  Et  [mt+i]  —  l/a  ■ 
1 

mib 


^     (mt  —  -),  so  that,  for  mt  in  a  neighborhood  of  1/a, 


Et 


1 

mt+i 

a 


p  \mt j  +oimt 

p    =     a/feG(0, 1). 


Hence  (34)  behaves  like  a  autoregressive  process  AR(1),  with  central  value  m,  =  1/a,  and  autocorre- 
lation coefficient  p  —  a/h.  However,  its  functional  form  yields  tractable  bond  and  perpetuity  prices, 
as  shown  by  the  following  Proposition.^^ 

Proposition  8    (Price  of  a  zero-coupon  bond  in  the  1-factor,  discrete  time  case)  The  price  of  a  zero 


There  are  many  ways  to  ensure  that  mt  remains  above  a  lower  bound  m  S  [1/6, 1/a).   The  simplest  one  is  to  get 
series  of  independent,  positive  random  variables  rj^,  with  E  [tjJ  =  1,  and  set;  mt+i  =  [^  +  j  —  ■^^ m)  rj  ^  +  rn.  As 

mt  >  me  [1/6,  1/a),  ^  +  ^  —  ^i^ m  >  0  and  rrtt+i  >  m.  The  reason  to  impose  mt  >  1/6  is  that,  if  we  had  mt  <  1/6 

for  some  t,  in  the  deterministic  version  of  the  process,  there  would  be  a  time  s  >  t  such  that  m.,  <  0,  which  is  not 
allowed  for  a  discount  factor. 

"The    correspondence    with    the    one-factor    process    is    the    following;      for    a    small    t,    a    =     1/(1— ra),  6    = 
1/  (1  -  T/3)  ,mt  =  1  -  rrt.  Then,  (34)  reads; 

rEt  [rt+i  -rt]  =  -  {E  [mt+i]  -  mi)  =  I mt)  {  -  -  mt]  —  =  (ri  -  q)  (n  -  /3) 

\a  J   \o  J  m.1  1  —  TTi 

i.e.,  E  [n+i  -  Tt]  =  {n  -a){n-l3)T  +  0  (r),  i.e.,  E  [dn]  =  (n  -  a)  (n  -  /3)  dt. 
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coupon  bond  of  maturity  T  is: 


_j^u  -  iilf  ,-T''^t       ~  "^ 


Zt  (T)  =  a-'  — '—  +  b-'  -i .  (35) 

0  —  a  b  —  a 


Summing  over  the  maturities,  we  get  the  price  of  a  perpetuity. 
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Proposition  9   (Price  of  a  perpetuity  in  the  1-factor,  discrete  time  case)  The  price  of  a  perpetuity 

is: 

„        a  +  b  —  I  —  m7  .„. 

^-=     (a-l)(6-l)-  •  (2«) 

The  above  discount  factor  can  be  enriched  by  a  stochastic  component,  as  in  Eq.  15. 

6.2     Discrete  time,  n  factors 

This  section  studies  the  case  with  a  number  of  factors  n  >  1.  A  quadratic  process  is  given  by: 
(i)  a  stochastic  discount  factor  mj  S  R*^_;  (ii)  a  state  vector  Xt  £  M";  (iii)  a  matrix  F  £  E"^", 
with  eigenvalues  in  (  —  1,1);  (iv)  a  vector  /3  €  R";  and,  (v)  a  central  value  of  the  gross  interest  rate 
i?*  =  1  +  r,  >  0,  with: 


Rtirit 


mt     =     ^{l-p'Xt). 

We  assume  that  the  process  is  well-defined  for  all  ^'s.  The  conditions  that  ensure  this  are  analo- 
gous to  the  continuous  case,  and  will  be  specified  in  the  next  iteration  of  the  paper. 

To  get  some  intuition  for  the  process,  we  first  consider  the  steady-state.  As  T  is  contractive,  the 
steady  state  value  of  X  is  0,  so  the  steady-state  value  of  mt  is  1//J,.  In  general,  if  T  has  positive 
roots,  then  Xt  is  positively  autocorrelated,  with  a  small  modulation  introduced  by  the  mt  term. 

Lemma  3  Suppose  that  process  yt  £  R'^satisfies  Et  [yt+i]  =  Cyt/mt,  for  some  qxq  matrix  C .   Then, 
vr  >  0, 

Et  [mt+i...mt+T-i2/t+r]  =  C^ — ■  (37) 

mt 

Proof.__Call  Lt-  the  left-hand  _side._I/i  =  _E\yt+\\  =.Cytlmt,  and 

Lt+i     =    Et  [mt+i...mt+Tyt+T+i]  =  Et  {mt+i...mt+T-iEt+T  [mt+TVt+T+i]] 
=    Et  {mt+i...mt+T-iCyt+T]  =  CEt  [mt+i...mt+T-iyt+T]  =  CLt- 


'   Alternatively,  one  can  try  a  price  process  of  the  type  Pi  =  A-\-  B  /mi,  and  solve  for  A,  B  that  satisfy  the  functional 
equation:  0  =  Pt  —  Et  [mt+i  (1  +  Pt+i)]- 
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Theorem  4  (Price  of  a  zero-coupon  bond  in  the  multifactor,  discrete  time  case)  The  price  of  a 
zero-coupon  bond  of ■  maturity  T  is,  with  /„  the  identity  matrix  of  dimension  n: 

Zt  (T)  =  r:^  -  R-J-'p'T  (/„  -  r)-i  (/„  -  r^)  ^.  (38) 

mt 

Proposition  10  (Price  of  perpetuity  in  the  multifactor,  discrete  time  case)  Assume  i?,  >  1.  The 
price  of  a  perpetuity,  Pt  =  '}Zt'=i  ^t  (■^)'  ^'''  ''^''■'l-^  ^n  the  identity  matrix  of  dimension  n: 

Proof.  One  simply  sums  (38).  Alternatively,  one  can  try  a  price  process  of  the  type:  Pt  = 
a  +  b'Xt/mt,  and  solve  for  a  6  K  and  6  e  R"  that  satisfy  the  functional  equation:  0  =  Pj  — 
Et[mt+iil  +  Pt+i)].    m 

7     Conclusion 

Quadratic  processes  are  very  tractable,  as  they  yield  closed  forms  for  perpetuities,  and  prices  that 
are  linear  in  factors.  They  are  likely  to  be  useful  in  several  parts  of  economics,  when  trend  growth 
rates,  or  risk  premia,  are  time-varying.  Hence  quadratic  processes  might  be  a  useful  addition  to  the 
economists'  toolbox. 
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Appendix  A.  Regularity  conditions  for  the  one-factor  process 

This  appendix  details  conditions  for  the  existence  and  uniqueness  of  the  solutions.  We  recommend 
Katatzas  and  Shreve  (1991  Chapter  5.5)  and  Revuz  and  Yor  (1999,  Chapter  IX)  for  systematic  treat- 
ments, and  Ait-Sahalia  (1996,  Appendix)  and  Dufhe  (2001,  Appendix  E)  for  pedagogical  overviews. 
We  call  V  =  (r,  r)  the  domain  of  existence  of  r,  and  c  an  arbitrary  point  in  V.  We  call  fi  (r)  the  drift 
of  r,  and  assume  dMt  =  a  (rj)  dzf.  We  make  the  following  assumptions. 

(i)  The  drift  and  diffusion  functions  are  continuously  differentiable  in  r  in  V,  and  a^  (r)  >  0  in 
D. 

(ii)  The  integral  of  m  (r)  =  exp  ( JJ^  2/i  {u)  ja"^  (u)  duj  /a'^  (r)  converges  at  both  boundaries  of  V. 

(iii)  The  integral  of  s  (r)  =  exp  (-  /J  2^  (u)  /a^  (u)  du)  diverges  at  both  boundaries  of  V. 

(iv)  ij,  is  Lipschitz  continuous,  and  there  is  a  function  p  {x)  :  R-|-  — >  R-j.,  with  p  (0)  =  0,  such  that 
for  any  e  >  0,  /,q  .  p {x)~   dx  =  +oo,  and  |cr  (x)  -  a (y)\  <  p{\x  -  y\) . 

If  conditions  (i)-(iv)  are  satisfied,  then  there  is  a  unique  Ito  process  {rt,t  >  0}  which  is  a  strong 
solution  of  the  stochastic  differential  equation  (1)  with  initial  condition  ro  =  r.  Moreover,  {rt,t  >  0} 
is  Markov. 

The  key  substantive  point  is  that  the  process  is  defined  for  all  i  >  0,  and  does  not  explode. 
This  condition  is  crucial,  as  if  we  started  with  ro  >  /?,  the  process  would  explode  in  finite  time  with 
positive  probability,  so  that  the  process  would  not  be  defined  for  all  times. 

Conditions  (i),  (ii)  and  (iv)  guarantee  the  existence  and  uniqueness  of  the  solution  up  to  the 
variable  may  hit  the  boundaries.  Condition  (iii)  implies  that  the  boundaries  are  actually  not  reached. 
The  intuition  is  as  follows.  Consider  the  correct  boundary.  Condition  (iii)  imphes  p.  (f)  <  0,  so  that 
the  process  tends  to  return  inside  P,  and  also  requires  that  a'^  (r)  tends  to  0  fast  enough  as  r  1  f . 

Sufficient  conditions  to  ensure  (i)-(iv)  Conditions  (i)  and  (ii)  guarantee  that  the  stochastic 
differential  equation  (1)  admits  a  unique  strong  solution.  Those  conditions  are  verified  in  the  following 
cases.  Condition  (iii)  guarantees  that  the  end  points  V  of  are  natural  boundaries. 

We  assume  p{f)  <  0  and  limr^r  p  (r)  >  0,  so  that  close  to  the  end  points  of  V,  the  process  tends 
to  go  back  inside  V.  In  the  case  p{r)  —  [r  —  a)[r  —  /3),  with  a  <  P,  this  corresponds  to  f  G  (q,/3) 
and  r  e  [— oo,  a). 

Conditions  (ii)  and  (iii)  are  verified  if  the  following  conditions  [C-V)  hold.  For  r  in  a  left- 
neighborhood  of  r,  cr^  (r)  ~  /c  (f  -  ry,  with  (k  >  1  and  fc  >  0)  or  (k  =  1  and  0  <  fc  <  -2m  (f)).  If 
r  >  — oo,  for  r  in  a  right  neighborhood  of  r,  a'^  (r)  ^  k'  {r  —  r)'^  ,  with  (k'  >  1  and  k'  >  0)  or  (k'  =  1 
and  0  <  fc'  <  2m  (r)).  If  r  =  —  oo,  then  r  is  a  natural  boundary  if,  for  r  in  a  neighborhood  of  —  oo 
_(j^(r)  ~  fc  \r\_,  w[th  k>  0  and  P  <  3,  Those  last  conditions  imply  assumptions  (ii),  (iiij.        


For  r  =  — oo,  the  situation  is  complex  for  condition  (iv),  as  the  standard  conditions  found  in 
textbooks  do  not  apply,  p  (r)  is  not  Lipschitz  continuous,  as  p'  [r)  is  unbounded.  We  conjecture 
that  a  simple  weakening  of  condition  (iv)  will  allow  the  case  r  =  —  oo. 

If  r  >  — oo,  the  above  conditions  (C-V)  also  imply  (iv),  as  one  can  take  p{x)  =  K  max  (x"'^,  x'^ '   ,x 
for  a  large  enough  constant  K. 
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Appendix  B.  Proofs 

7.1      Proof  of  Theorem  1 

It  is  sufficient  to  prove  the  statement  for  an  initial  time  equal  to  0.  We  fix  T  >  0  and  define 
Vt  =  Et    exp  (  —  Jj    r^ds)    ,  for  t  G  [0,T].  The  arbitrage  equation  of  Vt  is; 

-nVt  +  E  \dVt]  /dt  =  0  (40) 

We  seek  a  solution  Vt  =  Wt  with  the  functional  form: 

Wt  =  wirt,t)  =  A{t)rt  +  B{t)  (41) 

The  boiondary  condition  Wt  =  1  imphes  A  (T)  =  0  and  B  (T)  =  1 .  We  call  A  the  infinitessimal 
diffusion  operator  of  the  process,  which  in  the  case  dMt  =  cr  [rt)  dzt  is: 

Af  (r)  =  {r-a)ir-  P)  f  (r)  +  ^/"  (r) 

The  PDE  (40)  becomes: 

—rw  {r,  t)  +  Aw  (r,  t)  +  —w  {w,  t)  =  0 

Given  Ar  =  (r  —  a)  (r  —  /3),  we  have  Ati;  (r,  t)  =  A  (t)  (r  —  a)  {r  —  0).  Importantly,  because  d'^w/dr'^  = 
0  in  the  functional  form  (41),  the  Ito  second-order  term  (equal  to  ^—^d'^w/dr'^  in  the  case  of  a  dif- 
fusion) is  0.  The  PDE  is  Q  =  0  with: 

Q     =     —rw  (r,  t)  +  Aw  (r,  t)  +  —w  {w,  t) 

=     -r[A{t)r  +  B{t)]  +  A{t){r-a){r-p)+[A'{t)r  +  B'{t)] 
=     [aPA  (i)  +  B'  it)]  +  r  {-B  [t)  -  (a  +  /3)  A  (t)  +  A'  (t)] 

Importantly,  the  r^  terms  cancel  out.  This  explains  why  adding  a  quadratic  term  to  the  standard 
Ornstein-Uhlenbeck  process  substantially  improves  tractability.  The  PDF  Q  =  0  is  satisfied  for  all  r 
if  and  only  if: 

al3A{t)  +  B' {t)     =     0 
-^Bii)^-  (q  +-,0)-A(t)-+ -A^i)-  ^  -0 -(^)- 


Differentiating  the  last  equation,  we  get: 

-B'{t)-{a  +  p)A'{t)  +  A"{t)  =  Q 
which  imphes,  via  a^A  (t)  +  B'  [t)  =  0: 

a/3  A  (t)  -  (q  +  /3)  A'  (i)  +  A"  (i)  =  0, 
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which  is  satisfied  by  a  function  e^*  if  and  only  if: 

0  =  a/3  -  (a  +  /3)  e  +  ^'  =  (^  -  a)  (^  -  /3) 

so  that  the  roots  are  a  and  ^.    Therefore,  for  some  constants  c  and  C,  A  (i)  =  ce"'  +  Ce^'.    The 
equations  for  c  and  C  are  given  by  A  (T)  =  0,5  (T)  =  1  and  (42).  We  get: 


Wt  =  A  (t)  n  +  B  (t)  =  '- ^e"^'-^ '  + e 


P-a  P-a 


The  above  shows  that  Wt  satisfies  —rtWt  +  E[dWt]/dt  —  0,  and  Wt  =  1,  but  not  yet  that 
it  is  the  equal  to  the  bond  price.    To  conclude  the  proof,  we  define  Lj  =  exp  (  —  /q  Tsds\  Wt-    As 

Et  [dWt  -  nWtdt]  =  0,  we  have  Et  [dLt]  =  Et  [exp  (-  /J  r,ds)  {dWt  -  nWtdt) 

martingale,  which  implies,  Lq  =  Eq  [Lt],  i.e.  Wq  =  Eq   exp  ( —  j^  rgds  j  Wt 


=  En 


0.   Hence  Lt  is  a 


I.e. 


p  —  a  p  -  a 


exp 


r.ds 


which  is  the  statement  of  the  Theorem. 

7.2      Proof  of  Proposition  3 

The  Proposition  is  a  direct  re-expression  of  (3),  applied  to  the  process  rt  =  R  —  gt  which  follows, 
with  r^  =  R  -  5«:  drt  =  -4>{t*  -  r)  dt  +  {rt  —  r^)^  dt  —  dMt. 

We  also  present  an  alternative,  arguably  more  direct  derivation.  We  look  for  a  formula:  Vt  of  the 
type 

Vt  =  v  {Dt,gt)  =  Dt{a  +  b  {gt  -  g.)) 


which  implies: 


Dt 


E\dVt]/dt    =    E 


dDt 


Dt 


/dt  +  hE  \dgt]  /dt 
=    gt{a  +  b{gt- g,))  +  b  -4>{gt  - g*)  -  {gt  - g*) 
Expressing  the  above  in  terms  of  {gt  —  g*)  rather  than  gt,  the  PDE  is  Dt  —  RVt  +  E  [dVj]  /dt  =  0: 
_ --Q  -  - -J--  ^(a-+ 5 -(^j"^^))  +  [(^^  _  ^^y+  ^^]  [(a  +  i  (^-  -gM+  b  [(^  (gt  -  g»)  -  {gt  -'g*f 

=     [1-Ra  +  g,a]  +  {gt  -  5.)  [-Rb  +  a  +  bg,-  60]  +  {gt  -  g^f  [b  -  b] , 
which  is  satisfied  if  and  only  if: 

1  —  Ra  +  5,a    =    0 
-Rb-^a  +  bg^-bcf)    =    0. 

The  term  {gt  —  g*)    cancels  out,  which  justifies  its  inclusion  in  the  definition  of  the  quadratic  process. 
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The  last  two  equations  give: 


I.e. 


a    =     \/{R-g.) 

b    =    a/{R-g^  +  4>), 


v,  =  ,^(i+    ''-'* 


R-  g^  \       R-  g* 


7.3      Proof  of  Proposition  5 


This  is  a  simple  corollary  of  Theorem  1  and  Girsanov's  theorem.  To  simphfy  notation,  we  normalize 
the  initial  time  to  0,  g  =  0,  Dq  =  1.  Define: 

jj^  ^  eiv.-lN]./2 
where  IN]^  is  the  quadratic  variation  of  A''  up  to  time  t.  Define 

dM^  =  dMt-d  (Mt,  Nt)  =  dMt  -  [an  +  h)  dt 
then 

dn 


drt 


(r  -  A)  {r  -  id)  +  dMt 

{r'^  -  {X  +  IJ.-a)r  +  Xfi  +  b)dt  +  dM^ 

(r- A')  {r  -  n')  dt  +  dM^ . 


(43) 


Let  Q  denote  the  probability  with  Radon-Nikodym  derivative  dQ  =  e^'    '^"/•^dP: 


V  =  E 


^-J^'r^dujy 


-S^Tudu    Nt-\N]J2 


=  S<3 


-  /('  Tu  du 


By  Girsanov's  theorem,  dM^    is  a  martingale  under  Q,  so  Theorem  1,  applied  to  process  (43)  gives: 

7.4      Proof  of  Theorem  2 

It  is  sufficient  to  consider  X«  =  6  =  0,  by  redefining  Xt  to  be  Xt  —  X^  if  needed.  Also,  it  is  enough 
to  consider  r*  =  0,  by  redefining  rj  to  be  Tt  —  r,.  We  consider  the  Z  {Xt,  T)  the  price  of  zero-coupon 
bond  of  maturity  T.  The  arbitrage  equation  for  Z\s: 

0  =  -TtZ  +  E  [dZt]  /dt. 

We  guess  the  functional  form  Z  {Xt,T)  =  1  +  7'  (T)  Xt,  for  a  function  7  (T)  with  values  in  IR"  to 
be  determined.  With  this  functional  form,  E  [dZt/dt]  =  7'  (T)  {-AXt  +  nXt)  -  drZ  {Xt,  T),  so  that 
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the  PDE  becomes: 

0    =     -r,{\+^'{T)Xt)+i'{T){-AXt  +  rtXt)-^l{T)'Xt 
=     -n-j'mAXt-j^yiTyXt. 

The  rtXt  terms  cancel  out,  which  justifies  their  inclusion  in  the  quadratic  process.  Using  rj  =  (3'Xt, 
we  obtain: 

0     =     -/3'-y(T)A-^7(r)'  =  0 

^  jt^iT)'  =  -'y'(T)A-p'. 

One  integrates  this  matrix  equation  hke  a  scalar  equation,  using  the  boundary  condition  7  (0)  =  0, 
which  gives: 


7(T)'  =  /3' 


e-^^  -  1 


A 
and  finally  the  earlier  expression  for  the  bond  price. 

7.5  Proof  of  Proposition  7 

We  integrate  (22): 

Pt    =      /      Zt{T)dT=  e-'-'^dT+  e-'--^/?' (Xt-X,)dT 

Jo  Jo  Jo  A 

=     ^+P'  (f  e-('--+^)^  -  e-^'^dx)  J  (Xt  -  X.) 

=    -  +  P'(^r^--)U^t-X.)  =  --P'         \        (Xt-X.) 
r.  \r^+A       rtJA  r.  r^  [r^  +  A) 

7.6  Proof  of  Theorem  3 

We  use  the  convention  Jq  =  0  and  Sn+i  =  5n-  We  define  qj  =  1/(5,4. i  ^^'^ 


We  check  the  conditions  of  Lemma  2.  ip-    >  0,  and 


0  J  =  l 
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We  next  calculate  each  side  of  (30),  evaluated  at  the  k—th.  coordinate.  We  obtain: 


and,  using  J2j'^j    =  !> 


LHS,-1     =     (a.5+(l-a,^)rW)    -  1  = /^  +  fl  -  i^)  ^  -  1 
^i+i         J  \  Si 


RHS,-i   =   E^fr^^-E^;'  =  E^f  (ri'^-i 


J  3  3 


SiSi 


^J]<5,-(<5,+i-5,_i)l,<,(^-l 


If  k  >  L  RHSk  -1  =  0  =  LHSk  -l.llk<i, 


RHSk-1    =     yj—y]Sj{Sj+i-6j_,)l 


Sk      , 

k<j<i  \j:-i 


"^     k+l<i<i 


1 


5i5i+i 

1 

5i5i+i 
1 

5i5i+i 
LHSk  -  1 


\k+l<j<i 


k+l<j<i 


[(5j+i  +  5^-  5k+i  -  5k)  6k  -  6i+i5i  +  Sk+iSk 

1 


[((^i+i  +  5i  -  (5fc)  5k  -  5i+i5i 


5i5i+i 


{5i+i  -  5k)  {5k  -  5i) 


So,   we  get  V/c,  i^iJ^t  -  1  =  LHSk  -  1,  i.e.,  (30)  holds. 

7.7      Proof  of  Proposition  8 

We  have: 


Zt  (T  +  1)     =     Et  [mt+i...mt+T+2]  =  Et  [mt+i...mt+T+iEt+T+i  ['mt+T+2]] 

1 


Ef. 


1       1 

mt+i...mt+r+i  I  -  + 


a       b       abrnt-f-T+i 

-  +  -]  Et  [mt+i...Tnt+T+i] T-Et  [mt+i...mt+T] 

a       0  I  ab 


ZtiT  +  l)    =     {l  +  l)zdT)-^^Z,iT-l). 
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The  characteristic  equation  is:  ^^  =  (^  +  ^)  ^  -  ^,  which  has  roots  l/a  and  1/6.  Hence,  the  solution 
is  of  the  type:  Zt  (T)  =  Aa~^  +  Bb~^ ,  for  some  constants  A  and  B.  Using  the  values  Zt  (0)  =  1, 
and  Zt  (1)  =  Et  [mt+i]  =  l/a  +  1/6  —  1/  {abmt),  we  get  the  announced  formula. 

7.8      Proof  of  Theorem  2 

ZtiT)     =    Etlmt+i...mt+T]  =  Et[mt+i...mt+T-iK^{l-p'Xt+T)] 
=    R~^Et  \mt+i...mt+T-i]  -  K^P'Et  lmt+i...mt+T-iXt+T] 

r  V  Xt 


\R*J     mt 

as  the  previous  Lemma,   applied  to  Yt   =  Xf  and  C  =  F/i?,,  gives  Et  [Tnt+i...mt+T-iXt+T]   = 
(i^j'^^^-Thus, 


R'^Z  (T)  =  R^-^Zt  (T  -  1)  -  R^^pT^— 

mt 


and  so,  summing  from  1  to  T, 


RlZt  (T)  =  1  -  R-^P'  y  r^^  =  1  -  R-^/3'r^JL^^, 
■^       mt  In  -T  mt 

Appendix  C.  A  class  of  processes  admitting  a  closed  form  for  per- 
petuities 

The  following  Proposition  shows  a  way  to  generate  processes  that  have  closed  form  for  perpetuities. 
However,  typically  that  functional  process  does  not  yield  a  closed  form  for  finite-maturity  claims, 
unlike  quadratic  processes. 

Proposition  11  (Class  of  processes  admitting  a  closed  form  for  perpetuities)  Suppose  that  there  is 
a  process  iptt  <^''^d  constants  a  and  /3  such  that 

dipt  =  {rtipt  +  art  -  (3)  dt  +  dMt, 

where  dMt  is  an  adapted  martingale,   and  is  essentially  arbitrary  except  for  technical  conditions. 
'Then: 

V.  =  ^  (44) 

is  a  solution  of  the  perpetuity  PDE:  1  —  rtVt  +  E  [dVt]  /dt  =  0.  If  dMt  satisfies  regularity  conditions 
that  make  the  process  well-defined,  then  Vt  is  the  price  of  a  perpetuity,  Vt  =  Et 

Proof.  Call  Vt  =  (V't  +  a)  //3. 


J^oOe-/>udu^5 


1  -  rtVt  +  E  [dVt]  /dt=\-  rt'^^^j^  +  ^  (r^^t  +  o^n  -  /3)  =  0. 
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A  first  example  is  the,  quadratic  process:    Tp{r)  =  r,  a  =  —  (A  +  /x),5  =  —Xfi.    Eq.    44  gives 
Vt  =  (X  +  fi  —  rt)  /  (A/i),  the  formula  3  for  perpetuities. 
A  second  example  is 

d{l/rt)  =  <P{rt-r*)dt  +  dMt 

where  dMt  ensures  rt  >  0,  i.e.,  rt  mean-reverts  to  a  central  value  r*.  This  yields  a  price  for  perpe- 
tuities equal  to: 

F.=  (^  +  ^)/(l  +  ^)  (45) 

by  applying  Proposition  11  to  ip  (r)  =  1/r,  a  =  ^,b  =  l  +  (f)r* .  The  formula  (45)  admits  the  following 
interpretation.  If  rt  remains  at  its  current  value,  the  perpetuity  price  would  be  l/r-f .  If  rt  was  at  its 
long  run  central  value  r* ,  the  perpetuity  price  would  be  1/r*.  In  general,  the  perpetuity  price  (45) 
is  a  simple  weighted  sum  of  those  two  perpetuity  values,  (f)  indicates  the  speed  of  mean  reversion, 
and  gives  the  weight  between  the  "current  value"  1/r-j  and  the  "long  run  value"  1/r*. 
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